function Shell_clogit_FEdemoshrt_grID_wtp_tomlab_varagin_j_bs(subsample,inc,set_seed)

%matlab_seed=1353;
%rng(matlab_seed)

	tmp_file1=strcat('\\c3\rdat\SHoude\Research/sears/estar_data\refrigerators\choiceset_identifier_trimester_week_zipcode_2008_2012_v11022017_struct_',num2str(subsample),'_',num2str(inc),'_seed_',num2str(set_seed),'.csv');
    id=load(tmp_file1); 
	Nmax=size(id,1);
		if subsample == 11000
		bs_max=10*ceil(Nmax/10000);
    elseif subsample == 44000
		bs_max=10*ceil(Nmax/20000); 
    end

parfor j=1:bs_max	
	if subsample == 11000
		index_bs=[(1+(j-1)*(1000)):min(j*(1000),Nmax)];
    elseif subsample == 44000
		index_bs=[(1+(j-1)*(2000)):min(j*(2000),Nmax)];
	end
	Master_hom_FEdemoshrt_grID_wtp_tomlab_varagin_j_bs(subsample, inc, set_seed, j, index_bs)
end

end